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Abstract 

The Casimir energy or stress due to modes in a D-dimensional volume sub- 
ject to TM (mixed) boundary conditions on a bounding spherical surface is 
calculated. Both interior and exterior modes are included. Together with 
earlier results found for scalar modes (TE modes), this gives the Casimir 
effect for fluctuating "electromagnetic" (vector) fields inside and outside a 
spherical shell. Known results for three dimensions, first found by Boyer, are 
reproduced. Qualitatively, the results for TM modes are similar to those for 
scalar modes: Poles occur in the stress at positive even dimensions, and cusps 
(logarithmic singularities) occur for integer dimensions D < 1. Particular 

attention is given the interesting case of D = 2. 
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I. INTRODUCTION 



The dependence of physical quantities on the number of dimensions is of considerable 
interest ^$M. In particular, by expanding in the number of dimensions one can obtain 
nonperturbative information about the coupling constant |^|| . Useful expansions have also 
been obtained in inverse powers of the dimension J/J. 

In a previous paper we investigated the dimensional dependence of the Casimir stress on 
a spherical shell of radius a in D space dimensions [0. Specifically, we studied the Casimir 
stress (the stress on the sphere is equal to the Casimir force per unit area multiplied by 
the area of the sphere) that is due to quantum fluctuations of a free massless scalar field 
satisfying Dirichlet boundary conditions on the shell. That is, the Green's functions satisfy 
the boundary conditions 

G(x,t;x',t') =0. (1.1) 

|x|=a 

Here, following the suggestion of 0, we calculate the TM modes (for which H r = 0) in the 
same situation. The TM modes are modes which satisfy mixed boundary conditions on the 
surface [|T0|,|TT|j, 

|-r^ 2 G(x,t;x',t') =0, (1.2) 

Or |x[=r=o 

as opposed to the TE modes (for which E r = 0), which satisfy Dirichlet boundary conditions 
C P- . ID on the surface, and are equivalent to the scalar modes found in [§]. 

The organization of this paper is straightforward. In Sec. |T| we construct the Green's 
functions in Z)-dimensional space by direct solution of the differential equation, subject to 
the boundary condition (|1.2j ). Then, the Casimir stress is computed from the vacuum expec- 
tation value of the energy-momentum-stress tensor, expressed as derivatives of the Green's 
functions. The resulting expression for the Casimir stress on a D-dimensional spherical shell 
takes the form of an infinite sum of integrals over modified Bessel functions; the dimension D 
appears explicitly as well as in the orders of the Bessel functions. By combining the results 
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for the TM modes found here with those for the TE modes found earlier ||, we obtain a 
general expression for vector modes subject to perfect conductor boundary conditions on 
the spherical shell, which expression agrees with that found long ago for three dimensions 
fl"2|| . As a check, in Sec. [TT1] we rederive the same result from the vacuum expectation value 
of the energy density of the field. In Sec. [TV] we examine this expression for the TM Casimir 
stress in detail. We show that for D < 1 a constant can be added to the series without 
effect; a suitable constant is chosen so that each term in the series exists (each of the in- 
tegrals converges). We show how to evaluate the sum of the series numerically for all real 
D, using two methods: one involving Riemann zeta functions, and the second involving 
continuation in dimension. Both methods give the same numerical results. When D > 2 
the TM Casimir stress is real, and finite except when D is an even integer. The well-known 



D = 3 result |10|,|T3],^4|,|12j is reproduced when the n = mode is removed. When D < 2 
the Casimir stress is complex; there are logarithmic singularities in the complex- .D plane at 
D — 2, 1, 0, — 1, — 2, .... In the Appendix, the important case of D = 2 is discussed 

II. STRESS TENSOR FORMALISM 

The calculation given in this paper for the Casimir stress on a spherical shell follows very 
closely the Green's function technique described in ||, and we will therefore be concise, and 
emphasize the significant differences. 

A. Green's function 

The two-point Green's function G(x, t; x', t') satisfies the inhomogeneous Klein-Gordon 
equation 

(j^ - V 2 ) G(x, t- x', = -6^ (x - x')5(t - t'), (2.1) 
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subject to the boundary condition ( |1.2|) on the surface |x| = a. We solve this equation by 
the standard discontinuity method. In particular, we divide space into two regions, region I, 
the interior of a sphere of radius a and region II, the exterior of the sphere. In addition, in 
region I we will require that G be finite at the origin x = and in region II we will require 
that G satisfy outgoing-wave boundary conditions at |x| = oo. 
We begin by taking the time Fourier transform of G: 

/oo 
^e-^ ( *-*' ) G(x,t;x',t')- (2-2) 
-oo 

The differential equation satisfied by Q w is 

{uj 2 + V 2 ) &(x; x') = <^(x - x'). (2.3) 

To solve this equation we introduce polar coordinates and seek a solution that has cylin- 
drical symmetry; i.e., we seek a solution that is a function only of the three variables r = |x|, 
r' = |x'|, and 8, the angle between x and x' so that x • x' = rr' cos 9. In terms of these polar 
variables (EO) becomes 



a , & D-ld , sin 2 -^ d . D . 2Q d\ n , fCr-OWW) 



dr 2 r dr r 2 9^ 90/ v ' ' 7 2 7r( D - 1 )/ 2 r D - 1 sin D - 2 0' 

(2.4) 

We solve ( p.4|) using the method of separation of variables. The angular dependence is 
j;iven in terms of the ultraspherical (Gegenbauer) polynomial |L8| 



Ct 1+D/2 \z) (n = 0, 1,2,3, ...), (2.5) 



where z = cos 6. The general solution to ( |2.4| ) is an arbitrary linear combination of separated- 
variable solutions; in region I the Green's function has the form (with v = n — 1 + D/2 and 

k = M) 

oo 

g u (ry,0) =J2 a nr 1 - D/2 J»(kr)Ct 1+D / 2 \z) (r < r' < a) (2.6a) 

n=0 

and 



gUr,r',6) = Y,r l - D/2 [b n J v {kr) + c n .l. v {kr)\ C^ l+D ' 2 \z) (r 1 < r < a). (2.6b) 



n=0 



[Note that J v (x) and J_ v (x) are linearly independent so long as v is not an integer. Thus, in 



writing (|2.6b| ), we assume explicitly that D is not an even integer. We also assumed D > 2 



in writing down fl2.6a|) , so that J_„ is excluded because it is singular at r = 0.] The general 



solution to ( |2.4j) in region II has the form 

oo 

g^r,r ! ) 6) = Y,d n r 1 - D/2 H£\kr)Cl l + D/2 \z) (r > r' > a) (2.7a) 



n=0 

and 



g u {ry,0) = J2r X - D/2 [enH^\kr) + f n H<?\kr)] Cj- 1+D/2) (z) (r' > r > a). (2.7b) 

n=0 

The arbitrary coefficients a n , b n , c n , d n , e n , and f n are uniquely determined by six 
conditions; namely, the mixed boundary condition ( |1.2[ ) at r = a, 



D 
~2 



lj [b n J v (ka) + c n J- u (ka)] + ka[b n J' u (ka) + c n J'_ u (ka)] = (2.8a) 



^ - 1 ) \e n H^\ka) + f n Hl 2 \ka)} + ka[e n H^'{ka) + f n H®'(ka)] = 0, (2.8b) 



and 
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the condition of continuity at r = r', 

a n Ju(kr') = b n J u (kr') + c n J_ v (kr') (2.8c) 

and 

d n H£\kr>) = e n H^(kr') + f n H®{kr'), (2.8d) 

and the jump condition in the first derivative of the Green's function at r = r', 

b n J' v {kr') + c n J'_ v (kr') - a n J' u (kr>) = \| ; (2.8e) 

4{7rr') 2 k 

and 
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e n HW'{kr f ) + fnH^'(kr') - d n H^'(kr') 



(2.8f) 



4(irr') 2 k 

Solving these equations for the coefficients, we easily find the Green's function to be, in 
region I, 

g u {r,T>,9) = I ^(cosfl) [Mk r< )J^(k r> ) - (3J„(kr)J„(kr')} 



n=0 



(nrr') 



sm tcis 



(2.9a) 



where 



and, in Region II, 

G u (r,r',9) 



P 



(f - l) J- V {ka) + kaJ'_ u (ka) 
(f - 1) J v {ka) + kaJ' u (ka) 



E 

n=0 



2^ (f - 1) 

lG^rr')^ 2 " 1 ^ 



COS I 



H^\k r< )H^\h 



^H^ik^H^ikr') 



(2.9b) 



(2.10a) 



where 



(%-l)Hi 2) (ka)+kaH}, 2) '(ka) 



(2)// 
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(§ - l) Hi l \ka) + kaH ( y ] \k 



0)', 



(2.10b) 



B. Stress Tensor 



For a scalar field, we can calculate the induced force per unit area on the sphere from 
the stress-energy tensor T^(x, t), defined by 



T^(x,t) = <9V(x,t)<9<V(x,t) - l ~g^d x ip{x,t)d\{x,t). 



(2.11) 



The radial scalar Casimir force per unit area / on the sphere is obtained from the radial- 



radial component of the vacuum expectation value of the stress-energy tensor [12]: 



/=<o|ir-T-|o) u 



(2.12) 
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To calculate / we exploit the connection between the vacuum expectation value of the fields 
and the Green's function, 



(O|T0(x,t)0(x',t / )|O)=^(x,t;x / ,t / ) ; 



(2.13) 



so that the force density is given by the derivative of the Green's function at equal times, 
G(x,t;x',t): 



d d , , d d , . 

or or or or 



out 



(2.14) 



x=x', x =a 



It is a bit more subtle to calculate the force density for the TM modes. For a given 
frequency, we write 



(Trr) = - [V r V r , + U 2 - V± " Vj./] Qu. 



(2.15) 



where, if we average over all directions, we can integrate by parts on the transverse deriva- 
tives, 



n{n + D- 2) 

^2 ' 



(2.16) 



where the last replacement, involving the eigenvalue of the Gegenbauer polynomial, is ap- 
propriate for a given mode n [see (2.14) of ||]. As for the radial derivatives, they are[| 



V r = r z - u d T r 



n-2 



fl-io „AD-2 



(2-17) 



which, by virtue of ( |1.2| ), implies that the V r V r ' term does not contribute to the stress on 
the sphere. In this way, we easily find the following formula for the contribution to the force 
per unit area for interior modes, 



fi 



TM 
in 



7T 



{D+l)/2 2 D a D + l V (D^ 



r°° rln- J— , q 

^- / ^,(n,% 2 -n( n + D-2))^ 
—) Jo % n=(j s n 



X) 



(2.18) 



1 In the TM mode, the radial derivatives correspond to tangential components of E, which must 



vanish on the surface. See [11| 
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and for exterior modes, 
r . 

Jo 

where 



pTM 
'out 



n (D + l)/2 2 D a D+l T ( 



D-l- 



dx 



— ) J X n =0 



^2 w ( n i D )( x2 ~n(n + D- 2))— f 



w(n, D) 



(2n + D - 2)r(n + D - 2) 



n! 



(2.19) 



(2.20) 



x = ka, and the generalized Ricatti-Bessel functions are 



,( x ) = x^Ux), e n (x) = x D ' 2 - 1 H^\x). 



(2.21) 



It is a small check to observe that for D = 2 we recover the known result |T5 



eTM 
<D=2 



?7r 2 a 3 



oo 

E 



m 
x 2 



J m (z) , iZ^faO 



(2.22) 



where the half- weight at n = is a result of the /imi£ -D — > 2. In two dimensions, the vector 
Casimir effect consists of only the TM mode contribution. 

In general, we can combine the TE mode contribution, given in and the TM mode 
contribution, found here, into the following simple formula^]: 

w(n, D) 



f 



TM+TE 



E 



7r (r»+i)/2 2 r> a i?+i ^ r 

n=0 v 2 ' 

(s' n (x) e' n (x) <(x) 
Jo [s n {x) e n {x) s' n (x) e' n {x) 



(2.23) 



It will be noted that, for D = 3, this result agrees with that found for the usual electrody- 
namic Casimir force/area, when the n = mode is properly excluded. [See (4.7) of |I2| with 
the cutoff e = 0.] Of course, this only coincides with electrodynamics in three dimensions. 
The number of electrodynamic modes changes discontinuously with dimension, there being 
only one in D = 2, the TM mode, and none in D = 1, in general there being D — 1 modes. 



2 We will not concern ourselves with a constant term in the integrand, which we will deal with in 
Sec. IV. 
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Equation ( 2.23|) is of interest in a mathematical sense, because significant cancellations do 
occur between TE and TM modes in general. 



The integrals in (|2.18| ) and ( |2.19| ) are oscillatory and therefore very difficult to evaluate 
numerically. Thus, it is advantageous to perform a rotation of 90 degrees in the complex-a; 
plane. The resulting expression for / is 



./ 



TM 



n=0 



2 D ^a D 



w(n,D) [°° d 

/ ax x— in 

+i r (^i) Jo dx 



[2.24) 



III. ENERGY DERIVATION 

As a check of internal consistency, it would be reassuring to derive the same result by 
integrating the energy density due to the field fluctuations. The latter is computable from 
the vacuum expectation value of the stress tensor, which in turn is directly related to the 
Green's function, 

(Too) = ^ rfV + V-m (3-1) 

2 J-oo ^ r=r ' 

Again, because we are going to integrate this over all space, we can integrate by parts, 
replacing, in effect, 

V • V -> -V 2 -> oo 2 , (3.2) 

which uses the Green's function equation (|2.3[) . [Point splitting is always implicitly as- 
sumed, so that delta functions may be omitted.] Then, using the area of a unit sphere in D 
dimensions, 

2ir D / 2 

A ° = wiry (3 ' 3) 



we find the Casimir energy to be given by 

• 27r D/2 roc du 
T[Dj2)J-oo^ 



• 2 D/2 roc d ,oc 

E = ^ /n/0 , / -co 2 J r^drQ^r). (3.4) 



So, from the form for the Green's function given in ( |2.9a| ) and (|2.10a ), we see that we need 
to evaluate integrals such as 



r dr J v (kr)J- v (kr), 



which are given in terms of the indefinite integral 



(3.5) 



J dxxZ v {x)Z v {x) = y (Yl - ^ Z u {x)Z v {x) + Z' v {x)Z' v {x)\ , (3.6) 

valid for any two Bessel functions Z v , Z v of order v. Thus we find for the Casimir energy of 
the TM modes the formula 



E 



TM 



— — \ — — V w(n, D) 
2<kT{D - l)a t-i. 



dx 



x 



n=0 

2 ^+ D - 2 ))f5#l+^) + ( 2 -^ 



(3.7) 



x [ \^s^(x) e' n \x) 

We obtain the stress on the spherical shell by differentiating this expression with respect to 
a (which agrees with ( |2.18|) and ( |2.19| ), apart from the constant in the integrand), followed 
again doing the complex frequency rotation, which yields 



where the integrals are 



where 



q{x) 



oo 

f°° d 
Q n = - dxx—lnq(x), 
Jo dx 



y - l) Iu(x) + |(J v+1 (s) + Iu-xix)) 



(3.8) 



(3.9) 



x 



D 
~2 



x 



l)K v {x)--{K v+1 {x)+K v _ 1 {x)) 



(3.10) 



This agrees with the form found directly from the force density, ( |2.24| ), again, apart from a 
constant in the x integrand. 
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IV. NUMERICAL EVALUATION OF THE STRESS 



We now need to evaluate the formal expression ( |3.8|) for arbitrary dimension D. We 
implicitly assumed in its derivation that D > 2 and that D was not an even integer, but we 
will argue that ( |3.8|) can be continued to all D. 

A. Convergent reformulation of ( |3.8[ ) 

First of all, it is apparent that as it sits, the integral Q n in ( |3.9| ) does not exist. (The 
form in ( |2.24| ) does exist for the special case of D — 3.) As in the scalar case ||, we argue 
that since 



^2w(n,D) = for D < 1, (4.1) 

n=0 

we can add an arbitrary term, independent of n, to Q n in ( |3.8|) without effect as long 
as D < 1. In effect then, we can multiply the quantity in the logarithm in ( |3.9|) by an 
arbitrary power of x without changing the value for the force for D < 1. We choose that 
multiplicative factor to be — 2/x because then a simple asymptotic analysis shows that the 
integrals converge. Then, we analytically continue the resulting expression to all D. The 
constant —2 is, of course, without effect in fl3~9f), but allows us to integrate by parts, ignoring 
the boundary terms. The result of this process is that the expression for the Casimir force 
is still given by (|3.8|), but with Q n replaced by 



Q n = dx\n 







2 < \ 

■-q{x) 



(4.2) 



q(x) being given by ( |3.1U| ). 

Now the individual integrals in (|3.8| ) converge, but the sum still does not. We can see 
this by making the uniform asymptotic approximations for the Bessel functions in ( p,10|) 
fl9| , which leads to 

ixv ( -101 + 80D-16D 2 
Q n ~ — 1 + 



2 V 64i/ 2 

-5861 + 11152D - 768CLD 2 + 2304£> 3 - 256£> 4 

H 3 h • • • in oo). (4.3) 

16384^ 4 
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(Note that the coefficients in this expansion depend on the dimension D, unlike the scalar 
case, given in (3.17) of ||.) Because of this behavior, it is apparent that the series diverges 
for all positive D, except for D — 1, where the series truncates. 

There were actually two procedures which were used to turn the corresponding sum in 
the scalar case into a convergent series, and to extract numerical results, although only one 
of those procedures was described in the paper ||. In that procedure, we subtract from 
the summand the leading terms in the 1/n expansion, derived from (|4.3j ), identifying those 
summed subtractions with Riemann zeta functions: 
1 



F™ « 



1 N 

' n=l 



+ 



7T 



T(D-l) 



w(n, D)Q n — Tvn 

K+l 



N 



C(l - D) + h ^{k + 1-D)- b K+l 



n-K-2 



n 



(4.4) 



k=l n=l 

Here bk are the coefficients in the asymptotic expansion of the summand in (|3.8|) , of which 
the first two are 



6i 

^2 



(£>-2)(D-l) 
2 

81 - AA8D + 456D 2 - 17QD 3 + 2AD 4 
192 



(4.5a) 
(4.5b) 



In ( |4.4| ) we keep K terms in the asymptotic expansion, and, after N terms in the sum, we 
approximate the subtracted integrand by the next term in the large n expansion. The series 
converges for D < K + l, so more and more terms in the asymptotic expansion are required 
as D increases. 

There is a second method which gives identical results, and is, in fact, more convergent. 
The results given in || were, in fact, first computed by this procedure, which is based on 
analytic continuation in dimension. Here, we simply subtract from Q n the first two terms 
in the asymptotic expansion (f4.3|) , and then argue, as a generalization of (|4.1|) , that 



E 



T(n + D - 2) 



nl 



^T(n + D - 2) 9 
for D < 2, v . V 



for D < 0. 



nl 



(4.6) 



n=0 n=0 

Therefore, by continuing from negative dimension, we argue that we can make the subtrac- 
tion without introducing any additional terms. Thus, if we define 
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7T1S 



... , , -101 + 80D-16D 2 \ 

Qn = Qn ~ T ( 1+ ^ ) , (4.7) 



we have 



1 oo 

^™= 9 ' ^ w(n,D)Q n 
2na 2 T{D - 1) ^ 

I / N 00 

v 7 \n=0 n=JV+l 



r(n + D - 2) 
n! ^ 2 



(4.8) 



where ^(-D) is the coefficient of v 4 in (|4.3). The last sum in ( |4.8| ) can be evaluated according 
to 

y r(n + a) vr 2 r(a/2) 1 

^n!(n + «/2) 2 2 T(l - a/2) sin 2 vra/2' l ' J 



The approximation given in (|4.8|) converges for D < 4. 

B. Casimir stress for integer D < 1 

The case of integers < 1 is of special note, because, for those cases, the series truncates. 
For example, for D = only the n = 0, 2 terms appear, where the integrals cancel by virtue 
of the symmetry of Bessel functions, 

K v {x) = K_ u {x), I n {x) = J_ n (x), (4.10) 



for n an integer. However, using the first procedure (|4.4| ), we have a residual zeta function 
contribution: 

because both C(l — D) and T(D — 1) have simple poles, with residue —1, at D = 0. This 
result for D = is the negative of the result found in the scalar case, (3.22) of ||, which 
is a direct consequence of the fact that the n D ~ 2 term in the asymptotic expansion cancels 
when the TE and TM modes are combined [compare ( |4.5a| ) with the corresponding term in 



(3.23) of p.] The continuation in D method gives the same result, because then 
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F2i = + DQ, - 4) = ^ (x - ™ + «») = -L (4,2, 

where a limiting procedure, D — ► 0, is employed to deal with the singularity which occurs 
for n = 1, where i/ — > 0. 

For the negative even integers we achieve a similar cancellation between pairs of integers, 
with no zeta function residual because the ( functions no longer have poles there. For 
example, for D = —2 we have 

F™-2 = 7T^(Qo- 2Qi + 2Q 3 - Qi) = (4.13) 

because Q = and Q\ = Q 3 . Again, the other method of regularization gives the same 
result when a careful limit is taken. 

For odd integer < 1, trucation occurs without cancellation, because I„ ^ I_ u . For 
example, for D = 1, 

F ™i = 77-(<3o + Qi) = -0.2621 + 0.6032z. (4.14) 



C. Numerical results 

We have used both methods described above to extract numerical results for the stress 
on a sphere due to TM fluctuations in the interior and exterior. Results are plotted in Fig. 1. 
Salient features are the following: 

• As in the scalar case, poles occur for positive even dimension. 

• The integrals become complex for D < 2 because the function q(x), ( ^.10[ ), occurring 
in the logarithm develops zeros. (This phenomenon started at D = for the scalar 
case.) Correspondingly, there are logarithmic singularities, and cusps, occurring at 2, 
1, 0, —1, —2, . . . , rather than just at the nonpositive even integers. 

• The sign of the Casimir force changes dramatically with dimension. Here this is even 
more striking than in the scalar case, where the sign was constant between the poles 
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for D > 0. For the TM modes, the Casimir force vanishes for D = 2.60, being repulsive 
for 2 < D < 2.60 and attractive for 2.60 < D < 4. 



Also in Fig. 1, the results found here are compared with those found in the scalar or 
TE case, ||. The correspondence is quite remarkable. In particular, for D < 2 the 
qualitative structure of the curves are very similar when the scale of the dimensions in 
the TE case is reduced by a factor of 2; that is, the interval < D < 2 in the TE case 
corresponds to the interval 1 < D < 2 in the TM, — 2 < D < for TE corresponds to 
< D < 1 for TM, etc. 

Physically, the most interesting result is at D = 3. The TM mode calculated here has 
the value F™ 3 = —0.02204. However, if we wish to compare this to the electrodynamic 



result ||12|| , we must subtract off the n = mode, which is given in terms of the integral 
Qo = 0.411233 = 7r 2 /24, which displays the accuracy of our numerical integration. 
Similarly removing the n = mode (= — 7r/24) from the result quoted in (3.24) of ||, 
F™ 3 = 0.0028168, gives agreement with the familiar result ]10|,|13|,|IJ,|T2|,|] : 

0.0462 



pTM+TE 



9 ' 

d=3 a z 



(4.15) 



To conclude, this paper adds one more example to our collection of known results con- 
cerning the dimensional and boundary dependence of the Casimir effect. Unfortunately, we 
are no closer to understanding intuitively the sign of the phenomenon. 
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APPENDIX: TOWARD A FINITE D = 2 CASIMIR EFFECT 



The truly disturbing aspect of our results here and in [KJ are the pole in even dimensions. 
In particular many very interesting condensed matter systems are well-approximated by 
being two dimensional. Are we to conclude that the Casimir effect does not exist in two 
dimensions? 

One trivial way to extract a finite answer from our expressions, which have simple poles 
at D — 2 (I will set aside the logarithmic singularity there in the TM mode, because that 
only occurs in one integral, Qq), is to average over the singularity. If we do so for the scalar 
result in 0, we obtain 



while for the TM result here, we find 



0.340 



F^ 2 = (A2) 



which numbers, incidentally, are remarkably close to the leading Qq term, as stated^ in [15 



which are —0.0140 and —0.254. But, there seems to be no reason to have any belief in these 
numbers. 

However, something remarkable does happen in the scalar case. If we use the first 
procedure, fl4.4|) , we note that the poles can arise both from the integrals and from the 
explicit zeta functions. For the latter, let the dependence on D be given by r(D)/(D — 2) 
which has a pole at D = 2. When we average over the pole, we obtain 

lim i/ £ ^_ ! i2- £ n m 



-o2 

where the prime denotes differentiation. For the scalar modes it is easy to verify that 
r'(2) = 0. Thus, there is no contribution from those subtracted terms. In other words, they 



3 The integral in (A12) of [15| was not evaluated very accurately there. The value, good to 6 
figures, should be in our notation, -Q^ calar = 0.0880137. Similarly -Q™ = 1.5929. 
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might just as well be omitted, which is what we would do if we inserted a cutoff and simply 
dropped the divergent terms. (This procedure does give the correct D = 3 results.) This 
provides some evidence for the validity of the procedure which yields (|A1| ). 

Unfortunately, the same effect does not occur for the TM modes, r'(2) 7^ 0. Nor does it 
occur for higher dimensions, D = 4, 6, . . . , even for scalars. And, even for scalars, it is not 
clear how the divergences of a massive (2+1) theory can be removed. So we are no closer to 
solving the divergence problem in even dimensions.^ It is clear there is much more work to 
do on Casimir phenomena. 



4 For a discussion of the inadequacies of the dubious procedure of attempting to extract a finite 
result in [15 1 see Q. 
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FIG. 1. A plot of the TM Casimir stress F™ for -2 < D < 4 on a spherical shell, compared 
with F TE , taken from For D < 2 (D < 0) the stress F™ (F TE ) is complex and we have 
plotted ReF. 
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